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Abstract 

We investigate the emergence of Active Region NOAA 10790 by means of time- 
distance hehoseismology. ShaUow regions of increased sound speed at the location 
of increased magnetic activity are observed, with regions becoming deeper at the 
locations of sunspot pores. We also see a long-lasting region of decreased sound 
speed located underneath the region of the flux emergence, possibly relating to a 
temperature perturbation due to magnetic cjuenching of eddy difl^usivity, or to a 
dense flux tube. We detect and track an object in the subsurface layers of the Sun 
characterised by increased sound speed which could be related to emerging magnetic 
flux and thus obtain a provisional estimate of the speed of emergence of around 
1km s"^ 

Keywords: Sun: Active Region emergence; Sun: time-distance hehoseismology 



1. Introduction 

Local hehoseismology provides a tool for studying the emergence and evolution of 
active regions and complexes of solar activity. Several theoretical scenarios for the 
formation of active regions on the solar surface have been suggested within the 
context of the dynamo theory explaining the formation of the "sunspot belt" at 
low latitudes. Recently developed "flux-transport" dynamo models include the mag- 
netic flux transport by meridional circulation (Wang and Sheeley, 1991; Choudhuri, 
Schiissler, Dikpati, 1995; Dikpati and Charbonneau, 1999; Nandy and Choudhuri, 
2002; Krivodubskij , 2005), flux transport by an interface (Parker, 1993; Charbonneau 
and MacGregor, 1997) and overshoot dynamo with a positive radial shear beneath 
the surface (Ruediger and Brandenburg, 1995), or by a distributed dynamo with a 
negative shear (Roberts and Stix, 1972). 

All of the models agree that the most favourable and promising place for the aft 
dynamo is the region which includes the deepest layers of the convection zone, the 
convective overshoot layer and the tachocline. Nevertheless, it is difficult to keep ffux 
tubes with strong magnetic field (B >100 G) in the interior for a long time against 
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magnetic buoyancy (Parker, 1979) and then deliver tliem to the surface within the 
sunspot belt. 

In a flux-transport model it is possible to calculate the velocities of macroscopic 
diamagnetic transport as a function of depth in the convective zone by taking into 
account the magnetic advection processes. Krivodubskij (2005) estimated that in the 
subsurface layers the transport velocity increases from 0.5 km s~^ at a depth of 30 
Mm to around 0.9 km at 10 Mm depth. However, the way in which the flux 
travels to the surface is not fully understood. There are ongoing arguments in favour 
of each model, and they can only be tested by comparison witli the observations of 
active regions and sunspots on the surface and their appearance in the solar interior 
by means of helioseismology. 

Time-distance helioseismology (Duvall, Jefferies, Harvey et at, 1993) is a fast- 
growing area of helioseismology. The data used for time-distance inversions are the 
estimated wave travel-times between different points on the surface of the Sun. The 
travel-time estimates are obtained by cross-correlating and averaging the Doppler- 
velocity signal at different locations at the Sun's surface and then obtaining the cross- 
correlation function. The travel times are then computed from cross-correlations 
either via Gabor wavelet fltting (Kosovichev and Duvall, 1997) or by comparing 
with a reference cross-correlation function (see Gizon and Birch, 2005 and references 
therein). Because of the stochastic nature of solar oscillations, substantial spatial 
and temporal averaging of the data is required to measure the frequencies and travel 
times accurately. 

The properties of the solar interior are then determined by, flrst, establishing the 
relations between travel-time variations and internal properties (variations of the 
sound speed, flow velocity, magnetic fleld). This is followed by the inversion of these 
relations, which are typically cast as linear integral equations, to infer the internal 
properties. 

Chang, Chou, and Sun (1999) applied the metliod of acoustic imaging to study 
the emergence of NOAA 7978, and reported the detection of the signature of upward- 
moving magnetic flux in the solar interior. The flrst time-distance investigation of 
an emerging active region was by Kosovichev, Duvall, and Scherrer (2000), wherein 
an active region which emerged on the solar disk in January 1998, was studied with 
SOHO/MDI for eight days, both before and after its emergence at the surface. The 
results have shown a complicated structure of the emerging region in the interior, 
and suggest that emerging flux ropes travel very quickly through the depth range of 
our observations, estimating the speed of emergence at about 1.3 km s^^. 

In Section 2 of this paper the method and the data used in this work are described. 
In Section 3 the results of sound-speed inversions are presented, together with the 
discussion. This is followed by conclusions in Section 4. 

2. Data and Method 

2.1. Time-Distance Helioseismology 

The raw Doppler signal at any location is composed of velocity signals occurring due 
to convective motions on the scales of granulation and supergranulation, as well as 
the signal due to the superpositions of the wavepackets set up by resonant acoustic 
modes. The supergranular signal can largely be removed by flltering out frequencies 
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smaller than 1.5 mHz. However even after removing the low frequencies and frequen- 
cies above the acoustic cut-off of about 5.3 mHz, the oscillation signals due to the 
nature of sources and sinks of acoustic waves set up by turbulent convection remain 
stochastic. 

Due to the fact that acoustic waves with the same horizontal phase speed (to/k) 
travel the same horizontal distance (A) the Gaussian phase-speed filter is tuned to 
select waves with horizontal phase speed (vi) near the ray theory value corresponding 
to distance A^. We define a number of such filters (Fi) corresponding to different 
distances (A^) as follows: 

F,{k,u;;A,)=eM-{y-v,f/2Svf], 
k 

where cu is the frequency and k is the horizontal wavenumber in the data power- 
spectrum. Each filter is applied by pointwise multiplication of the Fourier transform 
S{kx,ky,u)) of the observed velocity data. 

Prom the filtered data, we calculate point-to-annulus cross-correlation functions, 
C{r, Ai,t), for each skip-distance and each position/pixel of our data, where r is the 
position vector, A^ is the distance between two surface bounces and t is time. The 
obtained set of cross-correlation functions is then averaged spatially over the quiet 
Stm to obtain an estimate for the reference cross-correlation function, C^°^(A,;,t). 
This is used to measure one-way travel-time perturbations at each location following 
Gizon and Birch (2002, 2004, also Jensen and Pijpers, 2003) approach: 

_ E.T/(±tOC'-°f(t.)[C(t.)-C"^(t.)] 

where f{t) is a window function designed to pick out the signal from the first bounce 
(i.e., direct arrival times). Here T-f and t_ correspond to wave travelling outwards 
and inwards from the measurement location respectively. We define the mean travel 
times r(r. A) as T(r, A) = (T+(r, A) + T_(r, A))/2. 

As a first approximation, perturbations [Srir, A)] in these mean travel times are 
linearly related to the sound-speed perturbations [5c] in the wave propagation region, 

Sr{r, A) = j j^dv' j ^ dzK{r - r', z; A)^(r', z), (2) 

where S is the area of the region, d is its depth and 5r(r, A) is defined as the 
difference between the measured travel time at a given location r on the solar surface 
and the average of the travel times in the cjuiet Sun. The sensitivity kernel for the 
relative squared sound-speed perturbation is given by K. 

For the forward problem we use wave-speed sensitivity kernels estimated using 
the Rytov approximation by (Jensen and Pijpers, 2003). The Rytov and ray approx- 
imations give very similar results, with the former being perhaps more reliable for 
inferring deeper structures (Couvidat, Birch, Kosovichev et al., 2004). 

The code for travel-time inversions is based on a multi-channel deconvolution 
algorithm (MOD) (Jensen, Jacobsen, and Christensen-Dalsgaard, 1998) enhanced 
by the addition of horizontal regularization. Following Jensen (2001) we discretized 
the Equation (2): di = (5T(k, A^), Gij = KCk, zj; Ai), rrij = Ss{]i,Zj). Then for 
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Table 1. Phase speed filter parameters 



i 




Skip-Distance range, Mm 


Vi, (km s '-) 


oVi, (km s j 


1 


5.83 


3.553-8.107 


14.21 


2.63 


2 


7.29 


5.013-9.567 


15.40 


2.63 


3 


11.664 


9.387-13.941 


17.49 


3.33 


4 


16.038 


13.761-18.315 


25.82 


3.86 


5 


20.412 


18.135-22.689 


30.46 


3.86 


6 


24.786 


22.509-27.063 


35.46 


5.25 


7 


29.16 


26.883-31.437 


39.71 


3.05 


8 


33.66 


31.383-35.937 


43.29 


3.15 


9 


37.91 


35.633-40.1870 


43.29 


3.15 


10 


42.28 


40.003-44.557 


47.67 


3.57 


11 


46.66 


44.383-48.937 


57.16 


3.78 


12 


51.01 


48.733-53.287 


57.16 


3.78 



each horizontal wavevector k we find the vector m that solves 

min{||(d-Gm)||l + e(k)2||Lm||i}, (3) 

where L is a regularisation operator. In this work we chose to apply more regularisa- 
tion at larger depths, where acoustic wave travel times provide less information, by 
setting L = diag(co,i, . . . , co.n), where co,j; is the sound speed in the i-th layer 
of the reference model. We heavily regularise small horizontal scales (Couvidat, 
Gizon, Birch et ai, 2005) by taking e(k)2 = el{l + (|k|/fcinax)^)^ with p = 200 
and eg = 5 X 10^. We invert for 14 layers in depth located at (0.36, 1.17, 2.11, 3.28, 
4.74, 6.41, 8.6, 11.22, 14.28, 17.78, 21.79, 26.31, 31.41, 36.95) Mm. 

The errors due to realisation noise are estimated by processing the 24 quiet- 
Sun datacubes, calculating noise covariance matrix in travel-time measurements 
(Jensen, 2001; Gizon and Birch, 2004; Couvidat, Gizon, Birch et al, 2005), and 
then evaluating noise in inversions according to (Jensen, 2001). 



2.2. Data 

We investigate the active region NOAA 10790 appearing on 13 July 2005. The full 
disk MDI data with one minute cadence are used. The time series start at 00:00 
UT, 10 July 2005, and ends at 23:59 UT 13 July 2005. From this set we extract 
forty four 512-minute series with the starting times every two hours. Each dataset is 
then re-mapped and de-rotated onto the heliographic grid using a Postel projection 
centred at 11° Carrington longitude and 11° latitude (South). The horizontal spread 
of the data is approximately 380 Mm. We use the Snodgrass rotation rate to remove 
the differential rotation (Snodgrass and Ulrich, 1990) for each processed datacube. 

Each of the datacubes is then processed to obtain the travel-time measurements 
using the centre-to-annulus cross-correlations as described in Section 2.1. The phase- 
speed filter values and the range of skip-distances used in this paper are provided 
in Table 1. The averaging annuli thickness was set to three pixels wide which 
corresponds to approximately 4.55 Mm. 
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Figure 1. The evolution of active region NOAA 10790 as observed by MDI. Carrington 
longitude is plotted along the x-axis, latitude along the y-axis. MDI continuum images are 
in the left column, the greyscale for MDI magnetograms (right column) corresponds to -600 
to 600 Gauss. The observation times for each row: a) 11 July 2005 11:11:32 UT, b) 12 July 
2005 06:23:32UT, c) 13 July 2005 06:23:32 UT, d) 13 July 2005 17:35:32 UT, c) 13 July 2005 
23:59:32 UT, f) 14 July 2005 11:11:32 UT. 
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Figure 2. Sunspot pore emergence. Top panels: the SOHO/MDI intensity (left) and mag- 
nctogram (right) images, tal^en on 12 July 2005 06:23 UT and 09:35 UT. Bottom panels: the 
cuts at 11 degrees of latitude through the sound-speed inversion estimated from 512 minute 
dopplergram timeseries centred on 11 July 2005 08:00 UT (left) and 11 July 2005 10:00 UT 
(right). The greyscale is in km s~^. Latitude is plotted on the left hand side along the j/-axis 
and longitude along the x-axis. Depth in Mm is marked on the right hand side of y-axis. 
Contours correspond to sound-speed perturbation values of 0.25 and 0.3 km s^^. 

3. Results and Discussion 

3.1. The Active Region NOAA10790 Surface Emergence and Evolution 

Signs of the increased magnetic activity at the site of the active-region emergence 
can be first observed on SOHO/MDI magnetogram images taken on 10 July 2005. At 
this time one can see the two points of flux at around 5° -8° Carrington longitude. 
The flux slowly grows witli the first pores appearing in the SOHO /MDI continuum 
images around 11 July 2005 06:30 UT (Figure 1). A number of different pores evolve 
and disappear at this location until 13 July 2005 06:30 UT, but no sunspots develop. 

Also around 14:30 UT 13 July 2005, we see a new flux emergence to the east of 
tlie first group at aroimd 14° longitude, 10-12° latitude (see Figure 1, panels d)- 
f)). At this new location, the first sunspot pores appear at about 17:30 UT 13 July 
2005, which quickly evolve into the bipolar sunspot group (by 14 July 2005 14:30 UT) 
before disappearing behind the limb. Also, tlrere is a bipolar region of relatively weak 
magnetic fiux seen appearing at around 20-22° longitude on 10 July 2005 and 12 
July 2005. No pores were observed at these locations. 

3.2. The Subsurface Sound-Speed Inversions 

Inversions of the sound-speed perturbation beneath tlie region of interest carried 
out in tliis work indicate that little change is detected at the first site of increased 
magnetic activity (6-10° longitude) prior to the emergence of magnetic flux at 
around 11:00 UT 10 July 2005. Then, once the magnetic flux appears on the surface, 
we see a consistent and substantial increase of the sound speed at shallow depths of 
up to 7-10 Mm directly beneath the locations of higli magnetic fleld strength, up 
until the end of the series used in this investigation. However, we have not been able 
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Figure 3. The sound-speed inversion of a 512 minute datacube centred around 13 July 2005 
04:00 UT. The top is the surface magnctogram at the start of the time series. The bottom is 
the surface magnctogram taken around the time of the end of the whole series (14 July 2005 
07:59 UT). The vertical panel is the absolute sound-speed perturbation at 11° latitude. The 
iso-surface corresponds to a value of 0.32 km s~^. 

to see any significant positive sound-speed perturbation in the deeper layers most of 
tlie time. Ttiis could be due to ttie fact tliat tlie spatial extent of tlie emerging flux 
at tfiis location is quite small, six -seven Mm in diameter for each footprint, and the 
signal cannot be resolved to better than half of the wavelength (Gizon and Birch, 
2004) . For the power-spectrum produced by the phase-speed fllters specified in Table 
1, this gives us an estimated resolution of, at best, 7.5 Mm for waves travelling to 
20 Mm depths. The averaging kernel at a target depth of 18 Mm, estimated at full 
width at half-maximum, extends from 8 to 22.5 Mm. This can also be a byproduct of 
the inversion regularisation, since due to our choice of operator, more regularisation 
is applied at deeper layers. 

As the ratio of gas pressure to magnetic pressure (the plasma /3) increases rapidly 
with depth, the sound speed increase due to magnetic field becomes very small, 
(5c/c of order 0.01 per cent at arotmd 18 Mm below the surface. Using the relations 
Pe = kTePe/m, pi = kTipi/m, = and pe = pi + pm, and assmning that the 
measured wave-propagation speed is given by = + c\, where the Alfven speed 
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Figure 4. The magnetic flux emergence. Tlie sound-speed inversion as a function of deptii 
and longitude. Tlie cut is made at 11° latitude. The inversion is obtained for the 512-minute 
datacubcs measured on 13 July 2005. The central times for the Dopplcrgram time series 
used to obtain the inversion are (top row): UT 04:00; 06:00; and (bottom row) UT: 08:00; 
10:00. Contours drawn for the features located between 4° and 16° longitude correspond to 
sound-speed perturbation values of 0.5, 0.7 and 1.0 per cent. 



is C4 = we can evaluate 

_ST 2_ ST 

Here p is pressure, T temperature, p density, k Boltzman constant, m mean molecular 
weight, and indices e,i,m stand correspond to external, internal, and magnetic; 
(3 = Pi/pm, &(? = — Ce, and 5T = Ti — Te. Thus, at greater depths the measured 
perturbation will be dominated by changes in temperature. This makes detecting 
and following localised perturbations due to emerging magnetic fields difficult. 

Notwithstanding the argument above, in our measurements the positive pertur- 
bation below the surface regions of magnetic activity becomes deeper and better 
defined as sunspot pores begin to form at the surface, indicating perhaps the increase 
of magnetic field strength. In the two bottom plots of Figure 2 we can see the region 
of the increased sound speed extending to depths of around 12-15 Mm at around 
7° longitude that, we believe, is related to the appearance of a sunspot pore seen 
in SOHO/MDI intensity images. By the time of 13 July 2005 00:00 UT there are 
several pores observed at this location, and the fiux configuration has a clear bi-polar 
structure (Figure Ic). 

Three-dimensional visualisation of the sound-speed inversion for an eight-hour 
series centered around 04:00 UT on 13 July 2005 is presented in Figure 3. One 
can clearly identify the subsurface tooth-like structure of the increased sound speed 
directly underneath the first region of surface fiux activity. In addition, an object 
characterised by increased sound-speed can be located deeper at the horizontal po- 
sition corresponding to the point of the second flux emergence. Longitudinal cuts 
in depth taken at 11° of latitude of sound-speed inversions for datacubes centered 
between 04:00 UT and 10:00 UT on 13 July 2005 are presented in Figure 4, with 
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Figure 5. Left: the results of the inversion at 21.77 Mm for the data cube around 13 July 
2005 04:00 UT. Right: results of the inversion at 6.41 Mm for the data cube centered around 13 
July 2005 08:00 UT. The contours in both plots are lines of ±150 Gauss from the magnetogram 
taken around 13 July 2005 17:35 UT. 



depth cuts at 21.8 Mm and 6.41 Mm of inversions for series centered on 04:00 UT 
and 08:00 UT plotted in Figure 5. 

One can see the region of increased sound speed at depths of 20 - 25 Mm at about 
14° longitude (see Figure 4, top left, and Figure 5, left). This region can be tracked 
rising on the next inversion as shown in Figure 4 (top right panel). The regions 
of increased sound speed can also be seen in the following two inversions (Figure 
4, bottom panels) rising to the surface at the location of the second flux group 
formation. The time of the appearance on the surface (13 July 2005 UT 06:00- 
14:30) corresponds to the first signs of the increased magnetic activity in the surface 
magnetograms. On the other hand, if the flux, by our assumption represented in the 
inversions as a region of the increased sound speed, travels from deeper layers to the 
surface, one should be able to see the traces of this flux presence at deeper layers in 
the inversions of eight-hour series centered at 08:00 and 10:00, while our inversions 
show only localised area of sound-speed increase. The averaged resolution kernels 
estimated for the depths layer around 21.79 Mm at half maximum extend from 8 to 
26 Mm. 

In order to check if the signals corresponding to sub-surface emerging flux add 
up while the noise structures cancel out, we have averaged all of our forty-four 3D 
sound-speed inversion maps. The averaging procedure is consistent within the linear 
limit used in the inversion scheme. The depth cut at 11° latitude of such an average 
is presented at the top of Figure 6, with daily averages plotted in other panels. 
One can clearly see the region of the increased sound speed corresponding to the 
flrst flux location. Interestingly, one can also observe a well deflned region of the 
decreased sound speed located beneath and to the West of the first fiux region with 
a column of negative perturbation corresponding to the location of the second fiux 
emergence. The structure is stable with respect to changes in the pre-processing 
(filtering) scheme and inversion parameters. It is also clearly present in the averages 
of the datacubes obtained on any of the days considered here, from 10 to 12 July 
2005 presented in Figure 6. However, the shape of this decreased sound speed cannot 
be observed in individual eight-hour inversions, possibly owing to turbulent activity. 
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We are not certain how to interpret this region of negative sound-speed perturba- 
tion but presumably it is a region of cooler plasma. Under certain conditions such a 
signal could arise from a monolithic flux tube originating in the deeper interior, with 
internal density exceeding the matter density outside the tube. Alternatively, such 
an effect can be seen in the scenario described by Kitchatinov and Mazur (2000) 
in which a temperature instability leads to the formation of an active region due 
to eddy diffusivity and magnetic tension. If the former explanation is the correct 
one, then it is tempting to identify such emerging cool flux tubes with the flux 
tubes postulated by Strous and Zwaan (1999) to explain transient darkenings they 
observed in continuum intensity. 

Thus, we can hypothesize that the sound-speed perturbation seen in Figure 4 is 
related to the emergence of a hot packet of plasma in the cold plasma column seen 
in Figure 6. In such circumstances, the absence of a positive perturbation column 
in the bottom two images of Figure 4 may be due to the dynamic nature of the 
region considered, wherein the body of increased sound speed preceded and followed 
by regions of decreased sound speed may cancel each other in our inversion. We 
note that the inversions of the interior sound speed averaged over eight datacubes 
obtained on 13 July 2005 (Figure 6, bottom right panel) point to the presence of a 
decreased sound-speed structure column at the location of the second flux emergence 
with a positive perturbation directly below the surface. The latter feature can be 
expected given that these data are for a period after the flux emergence at the surface 
has become well established. So, if we assume that the region of the increased sound 
speed corresponds to the emerging magnetic flux and compare the sound-speed per- 
turbation detected at the depth of 21.7 Mm on 13 July 2005 04:00 and appearing on 
the surface around 10:00 UT, this can provide an estimate of the speed of emergence 
at about 1.0 kms~^. This estimate is slightly lower than the speed reported by 
(Kosovichev, Duvall, and Scherrer, 2000), but deflnitely higher than the velocities of a 
diamagnetic transport in near-surface layers estimated in (Krivodubskij, 2005). This 
can be an indication of some additional mechanisms affecting the flux emergence. 

In our solutions, the errors due to realisation noise are less than 100 ms~^: this 
is estimated by processing the 24 quiet-Sun datacubes and applying the approach 
described in (Jensen, 2001; Gizon and Birch, 2004; Couvidat, Gizon, Birch et al., 
2005). 

4. Conclusions 

We have investigated the emergence of the Active Region NOAA 10790. The region 
has a complex magnetic structure and can be viewed as an amalgamation of two 
magnetic region as described in Section 3.1. We see shallow regions of increased 
sound speed at the location of increased magnetic activity, with regions becoming 
deeper at the locations of sunspot pores. It appears that a subsurface sound-speed 
perturbation at depth is associated with the emergence on the surface of a sunspot 
and sunspot pores. We observe an object with increased sound speed which may 
be related to the subsurface flux prior to its emergence in the photosphere at the 
second location. We also observe a loop-like structure of the decreased sound-speed 
perturbation located below the investigated region, which is present on the averages 
of the inversions. It would be interesting to see if similar structures characterised by 
a drop in sound speed can be detected in other instances of emerging active regions. 
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Figure 6. The mean sound-speed perturbation, averaged over all of the processed data, 44 
datacubes, (top left) and the mean sound-speed perturbation for the selected dates (averages 
of 12 datacubes). The greyscale is (5c/c in percent. The contours correspond to sound-speed 
perturbation of 0.1 percent. 



Based on our results, one can speculate that magnetised regions without any 
sunspots or sunspot pores lead to increased subsurface sound speed at shallow depths 
but no perturbation detected at greater depths. On the other hand, the areas associ- 
ated with sunspot appearance have increased finger-like sound-speed perttirbations 
which extend to depths of 10-20 Mm. Then by tracking the times of the emergence 
of magnetic flux in the subsurface layers of the Sun one can estimate the speed at 
which flux emerges to be about Ikms^^. 
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We also note that the precise timing of the signature of emergence in the seismic 
imaging is very hmited by the temporal resolution of the sound-speed inversions 
presented here, since each set is a representative of eight hours of solar oscillations. 

In view of the dynamo models discussed in Section 1 (Nandy and Choudhuri, 2002, 
Krivodubskij , 2005), it could be interesting to compare such results for a variety of 
emerging active regions at different latitudes and stages of the solar cycle. 

Acknowledgements We would like to thank an anonymous referee for the valu- 
able suggestions that have greatly improved this paper. We also acknowledge the sup- 
port from UK Science and Technology Facilities Council grant number PP/C502914/1.| 



References 

Chang, H.-K., Chou, D.-Y., Sun, M.-T.: 1999, Astrophys. J., 526, 53. 
Charbonneau, P., MacGregor, K.B.: 1997, Astrophys. J., 486, 502. 
Choudhuri, A.R., Schiissler, M., Dikpati, M.: 1995, Astron. Astrophys., 303, 29. 
Couvidat, S., Birch, A.C., Kosovichev, A.G., Zhao, J.: 2004, Astrophys. J., 607, 554. 
Couvidat, S., Gizon, L., Birch, A.C., Larsen, R.M., Kosovichev, A.G.: 2005, Astrophys. J., 
158, 217. 

Dikpati, M., Charbonneau, P.: 1999, Astrophys. J., 518, 508. 

Duvall, T.L. Jr., Jefferies, S.M., Harvey, J.W., Pommerantz, M.A.: 1993, Nature 362, 430. 
Duvall, T.L. Jr., Kosovichev, A. G., Scherrer, P. H., Bogart, R. S., Bush, R. I., de Forest, C, 

Hoeksema, J. T., Schou, J., Saba, J. L. R., Tarbell, T. D., Title, A. M., Wolfson, C. J. and 

Milford, P. N. : 1997, Solar Phys. 170, 63. 
Gizon, L., Birch, A.C.: 2002, Astrophys. J. 571, 966. 
Gizon, L., Birch, A.C.: 2004, Astrophys. J. 614, 472. 
Gizon, L., Birch, A.C.: 2005, Living Reviews in Solar Physics 2, 6-. 
Jensen, J.; 2001, Ph.D. thesis. Univ. Aarhus. 
Jensen, J., Pijpers, P.P.: 2003, Astron. Astrophys. L5, 461. 

Jensen, J., Jacobsen, B.H., Christcnscn-Dalsgaard, J.: 1998, In S. Korzcnnik & A. Wilson 

(eds.) Proc. SOHO 6/GONG98 Workshop, SP-418, ESA, Noordwijk, 635. 
Kitchatinov, L.L., Mazur, M.V.: 2000, Solar Phys., 191, 325. 

Kosovichev, A.G., Duvall, T.L. Jr.: 1997, In P.P. Pijpers, J. Christensen-Dalsgaard, and C.S. 

Rosenthal (eds) SCORe'96 : Solar Convection and Oscillations and their Relationship, 

ASSL 225, Kluwer Academic, 241. 
Kosovichev, A.G., Duvall, T.L. Jr., Scherrer, P.H.: 2000, Sol.Phys., 192, 159. 
Krivodubskij, V.N.: 2005, Astron. Nachr., 326, 61. 
Nandy, D., Choudhuri, A.R.: 2002, Sci., 296, 1671. 

Parker, E.N.: 1979, Cosmical Magnetic Fields, Clarendon Press, Oxford. 

Parker, E.N.: 1993, Astrophys. J., 351, 309. 

Roberts, P.H., Stix, M.: 1972, Astrophys. J., 18, 453. 

Ruediger, C, Brandenburg, A.: 1995, Astrophys. J., 296, 557. 

Snodgrass, H.B., Ulrich, R.K.: 1990, Astrophys. J., 408, 707. 

Strous, L.H., Zwaan, Z.: 1999, Astrophys. J., 527, 435. 

Wang, Y.-M., Sheeley, N.R. Jr: 1991, Astrophys. .J., 375, 761. 



sz_EAR10790.tex; 18/07/2008; 16:09; p. 12 



